Conditions for Factorizable Output Prom a Beam splitter 
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A beam splitter is one of the most important devices in an optics laboratory because of its 
handiness and versatility; equivalent devices are found in various quantum systems to couple two 
subsystems or to interfere them. While it is normal that two independent input fields are superposed 
at the beam splitter to give correlated outputs, identical Gaussian states interfere there to produce 
totally independent output fields. We prove that Gaussian states with same the variance are the 
only states which bring about factorizable output fields. 

PACS numbers: 42.50.-p, 42.50.Dv 



I. INTRODUCTION 

Quantum state preparation is the manipulation of a 
certain state into a desired target state. It is crucial 
for both the testing and implementation of quantum me- 
chanics to be able to generate custom states (or custom 
states within a specified tolerance) 

There are two main approaches to quantum state 
preparation: either by the time evolution of the state 
using a Hamiltonian that transforms the state from its 
initial condition to the final target state or by conditional 
preparation. In the second approach a measurement is 
made on one part of a bipartite correlated system; the 
action of this measurement is to modify the correlated 
system with the aim of producing a useful state. Paris et 
al. discuss the conditional preparation of a state in one 
mode of a two mode squeezed state by an appropriate 
measurement on the other mode in Ref. [2j • Conditional 
preparation is an essential feature of teleportation proto- 
cols HQ. 

An important application of the conditional prepa- 
ration method is the Knill-Laflamme-Milburn (KLM) 
scheme [5] for optical quantum computing based on the 
use of measurements to evolve states by exploiting sym- 
metry relations satisfied by bosons. This work was sem- 
inal as it theoretically demonstrated the possibility of 
scalable linear optical computing. Previous implemen- 
tations of quantum circuits had been proposed which 
were linear, Cerf et al. Q had put forward an all lin- 
ear scheme that avoided introducing the non-linear Kerr 
effects which were a common feature of schemes up until 
then (for example Q)- The KLM scheme does not suf- 
fer from the scalability problems of the other proposals 
which required an exponential increase in the number of 
optical elements when processing increasing numbers of 
qubits. Much work has been done since, to refine this 
protocol and improve its efficiency Q. 

For continuous-variable fields, it was proven Q that a 
quantum state can be engineered by single photon ad- 
ditions together with displacements [lfj. As the pho- 



ton addition is relatively difficult [Tl|, Fiurasek et at. 
[l2| came up with an idea to subtract photons from a 
squeezed field to generate a desired quantum state. Pho- 
ton subtraction has been realized in various laboratories 
using a beam splitter and a single photon detector (l3| . 

In quantum state engineering, two output fields are 
in general (quantum-mechanically or classically) corre- 
lated. In this way, by measuring one field, the other field 
is collapsed into a desired state. It was shown that when 
nothing is injected into one input port of a beam split- 
ter, the other input has to be a coherent state for the 
two output fields to be factorizable [14[. However, we 
can question if there are general input states which give 
totally uncorrelated outputs. In this paper, we seek to 
answer this question. 

We begin with a brief review of the topic of beam 
splitters (Sec. |TTJ) and the ability to perform condi- 
tional preparation with only classically correlated states 
(Sec. lII A"|) . In Sec. lIIII we examine the conditions required 
for the output of a beam splitter to be factorizable, and 
hence not correlated. The result that factorizable output 
requires Gaussian states as input is presented and this 
leads us to specify that, in fact, the Gaussian states may 
differ only by an arbitrary displacement. An alternative 
description, in the Wigner function representation, is also 
examined. 



II. STATE PREPARATION WITH A BEAM 
SPLITTER 

Beam splitters arc ubiquitous in quantum optics and 
are an essential component of almost all photonic quan- 
tum state engineering protocols. They are passive (en- 
ergy conserving) devices which are composed of a linear 
medium; that is the polarisation is proportional to the 
incoming field with constant of proportionality given by 
the first order susceptibility ■ They may be two at- 
tatched glass prisms or a partially silvered mirror and 
are constructed to a given specification so as to transmit 
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a certain proportion of an incoming source and reflect a 
corresponding amount as quantified by the coefficients of 
transmittivity and reflectivity (t — cos | and r = sin | 
respectively and satisfying r 2 + t 2 — 1). For two inci- 
dent beams the beam splitter facilitates interference of 
the two input fields to create two, generally different, 
output fields. 

Quantum mechanically we describe the action of a 
beam splitter by the operation [Til ] 



B(9) = exp[i-(tfb + atf)] 



(1) 



where d and b are the annihilation operators of the two 
input modes and d^ and w are their Hermitian conju- 
gates. The need for the device to preserve the bosonic 
commutation relations in the output modes prescribes 
that the operator must be unitary. The trivial cases of 
r = or t = are ignored since in these situations the 
beam splitter is not mixing the two input modes and so 
no correlations may arise. 

In the Heisenberg picture the action of the beam spitter 
is to effect the following changes: 



d — > ta + rb 



tb- 



(2) 



Defining quadrature operators q a — -j=((V + a), p a — 
-^(d^ — d), and similarly q b and p b , the relations describ- 
ing the first order transformations ([2]) may be rewritten 
in terms of the quadratures qi and pi as 



(3) 
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Using the definition of the two-mode Glauber-Sudarshan 
P-function (P function for simplicity) [13] 

P = J P(a,f3)(\a) aa (a\®\P) bb (f3\)d 2 ad 2 [3 (6) 

and Eq. (O we easily find the P function for the output 
field to be 
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(7) 



This P function is a smooth, non-negative function and 
thus the state §5§ describes a state which is classically 
correlated [2(| • Despite the lack of quantum correlations 
this classical state may be used in an implementation 
of conditional measurement. For example, if a measure- 
ment results in the second mode being projected onto the 
one photon Fock state then mode a becomes the non- 
Gaussian state 
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where ^4=1 
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surements on one mode of a classically correlated bipar- 
tite state a non-Gaussian state has been engineered in the 
other mode. Thus, even classically correlated states may 
be useful in quantum state engineering by the method 
of conditional measurement. Here we have produced a 
non-classical state 16] 



III. FACTORIZABLE OUTPUT FROM A BEAM 
SPLITTER 



A. State Modification by Measurement of a 
Classically Correlated State 

To implement conditional preparation requires two 
correlated states, however this correlation need not 
be entanglement — classical correlations can be suffcient. 
Consider the case of two distinct thermal states which 
represent equilibrium states for given temperatures Ti 
(i = a,b): 
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(4) 



where rij is the mean number of photons and \a) rep- 
resents a coherent state of amplitude a. Impinging two 
distinct thermal fields p a and p b onto a beam splitter 
generates the state with density operator 
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In this section, we find the condition for factorizable 
input fields of a beam splitter to result in completely 
independent output fields. Such states can be of no use 
in quantum state preparation. 

A density operator, p, for continuous-variable fields is 
assumed to be written as the exponential of a Hermitian 
operator / as 



P : 



J(a,al 



(9) 



Continuous- variable fields necessarily involve unbounded 
operators such as the quadratures q and p. Such un- 
bounded operators have spectra that may go to infinity 
and this fact implies that unbounded operators should 
be definied on certain subspaces of the Hilbert space due 
to the Hellinger-Toeplitz theorem [T3|. For these cases, 
every functional or density operator can be defined by 
employing unbounded operators such as d and d^ on their 
intersected subspace from spectral theory . The den- 
sity operator of Fock states and superpositions of Fock 
states can be written in the form (O . This is not the case 
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for all operators on a finite dimensional Hilbert space: A 
finite-dimensional density operators may have zero eigen- 
value (s) and so may not be written in such an exponential 
form with the exponent being any bounded operator; ex- 
amples would be density operators corresponding to spin 
states, polarisation states and angular momentum states. 

Equally we may define the functionals p and / in terms 
of the quadratures q and p as 



(10) 



We shall represent the exponent / by a power series with 
respect to q and p. The expansion could be made with 
respect to any ordering. As the exponent / is Hermi- 
tian, we choose symmetric ordering, denoted by {-} s ; 
for example, {£i£j} a = ^(CiCj + and {i, i', 6 |> = 
^7 (sum of all permutations of ^^j^fc), where £ = (q,p). 
Then, 

/ = /(o) + + &•}. + f${£i £j in}. + ■■■ 

(ii) 

where i,j, k, . . . — 1,2. Here we employ Einstein's con- 
vention for the summations. Because of the ordering 
used, the n-th order coefficient tensor /(„) is of all per- 
mutations consisting of n-th moments in q and p. The 
fact that the exponent / is Hermitian and the expansion 
is made in the symmetric ordering in terms of q and p 
implies that all the coefficient tensors /(„) are real and 
symmetric. The coefficient tensors /(„) completely spec- 
ify the state. The representation in Eq. (Ti"Tj) contains 
redundancies in the sense that all components of /(„) 
are not independent as /(„) is symmetric; for example, 
= ffk ■ Remembering this redundance, we shall still 
employ the expansion in the form of Eq. pip for mathe- 
matical convenience. 

The transformation rules for f( n \ follow once those for 

£j are known. Supposing that £j is transformed by ap- 
plying a unitary operation as 



(12) 



we then have the transformation rules for the (symmet- 
ric) higher moments, for example, 

u&ijikhfr = AlAfAMduin}.. (13) 

The transformed operators, U^iU\ preserve the commu- 
tation relations of £j. The transformed density operator 
obtained by applying the unitary operation is now rep- 
resented as 

Upfr = exp(#/f> t ). (14) 
Here, the transformed exponent is expanded as 

UfU ] = /(o) + /]d £ + /|) {£i£jh + ---, (15) 



where the transformed coefficient tensors /(„) are given 
in terms of /(„) and A^ as 
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Thus the action of an arbitrary unitary operation pre- 
serves the form of p — e^'^ by only evolving the coeffi- 
cients of / (coefficients get multiplied by the A^ effecting 
the transformation of the symmetric products which in 
turn depend only on the first order/linear transforma- 
tions). 

Now, consider a factorizable two-mode state, 



p apb = e /«.rf.)+9(ftA) < 



(19) 



The exponent of the two-mode density operator p a b — 
p a Pb is given by h = / + g and can be expanded in terms 
of £= {q a ,Pa,qb,Pb) as 



h = h {Q) + h l (1) ii + h^ 2) {£i £j} s + 



(20) 



where the coefficients h results from those of / and g, 
that is, 

h(o) = f(o) + 9(0) 
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It is clear that the factorizable state is conditioned by 
the vanishing coefficients of the mixed moments such as 
/i(2) and h l fo\ . Note that the mixed moments appear from 
second order onward. For the moments higher than first 
order, all the coefficients h % ^?"' %n are of mixed moments 
except the cases where Vife < 2 and Vz/c > 2. 

The transformations by the beam splitter of the co- 
efficient tensors hi n ) are given by Eq. ([3]), similarly to 
Eqs. (TT2]) - (HHJ) - The first order coefficients are trans- 
formed to 



h l — A* h 3 

n (l) ~ 1V 3 n (l) 



where 
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In general, the transformations of the n-th order coeffi- 
cients hr n ) are given as 



(n) 



A' 1 Af 



A' 



'/iff 

, (n) 



(24) 



In order for output to be factorizable/uncorrelated, all 
the vanishing coefficients in Eq. (|2Tj) must be retained 
after the transformation. For factorizable output we have 
no restriction on the zero-th and first order moments. For 
the second order moments, we have the conditions, 



p j+2 



0. 



(25) 



for each i, j < 2. The other conditions h 1 ^, J = are re- 
dundant due to the symmetry of coefficient tensor h(2)- 
The conditions (|2"5|) are rewritten in terms of the trans- 
formation matrix A and the untransformed coefficients 
h(2) as 



AlAf 2 



hf 2) = 0. 
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Similarly, for j < 2, A:J +2 = r 
if k ^ j or j 

or both k, I > 2. When these are applied, the conditions 
become 



2. Further h'fL ^ only if both k, I < 2 
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The case of r 
ther simply a 
the case of t 



= or t = is irrelevant as it implies ei- 
mirror or no beamsplitter, and we consider 



and J+2 



^ and r ^ 0. Noting that /i!-L = 



/(2) 



5^2) f° r hj ^ 2 (when the inputs are 
factorizable), the conditions imply that 



J{2) ~ 9(2) 



(28) 



This is the condition for the output to be factorizable for 
factorizable input states that are Gaussian (i.e. second 
order). As we shall now see in fact it is only second 
order input states that are permitted if the output is to 
be uncorrelated. 

Similar arguments can be applied to the higher order 
coefficient tensors The factorizable conditions, sim- 
ilar to Eq. (|2"T|) , are given as 

rh n-j f i^---i n + ( _ r) 



0. 



(29) 



for j = 1, 2, . . . ,n — 1 and all ik- Assuming t ^ and 
r 7^ 0, the conditions are reduced for odd n-th order 
moments into 
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which clearly result in, for all ik 
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Similarly, for even n-th order moments, the conditions in 
Eq. (f2"9"|) are reduced to 
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tions (|3Tj) hold for all moments higher than second order. 

All terms above second order must be zero in an in- 
put state in order for the output from a beam splitter to 
be factorizable. Furthermore the conditions in Eq. (|28|) 
specify that the two states should have the same second 
order coefficients whereas the first order terms are not 
restricted, i.e. Gaussian inputs differing only by a dis- 
placement give factorisablc output. 

Let us examine why they may differ by a displacement. 
Arbitrary displacements of n-th order input states creates 
additional terms of order n — 1 and lower which may then 
be subsumed into the coefficients of existing terms. Con- 
sidering the restriction to second order the displacements 
may only introduce linear terms and these terms cannot 
alter the factorisability of the output (these linear terms 
when transformed, according to Eq. ([3]), remain linear). 

Alternatively, we may see this condition in the Wigner 
representation. Based on the neccessity of cross terms 
canceling in a general form of the Wigner function, we 
can displace identical Gaussian input states each by arbi- 
trary amounts and still retain factorizability of the out- 
put. A single-mode Wigner function is given by [2l[ 

''/• 



W(q,p) 



2- 



Z l ™(q- V -\p\q+ V -)dy (33) 



where \q) are eigenstates of a quadrature variable q. The 
Wigner function for a general Gaussian state is the ex- 
ponential of a quadratic function of z = (q + ip) . One 

form, of Argawal [22j , that lends itself to describing many 
states is 

Wx(z) ' 



exp 
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which corresponds to the density matrix 
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(35) 

(36) 

(37) 
(38) 

(39) 



5 



The state is entirely specified by p, r and zq where z de- 
scribes the degree of displacement from the vacuum state. 
The requirement that a density matrix corresponds to a 
legitimate physical state (i.e. is positive semidefmite) is 
expressed as 



(40) 



in terms of the parameters r and p [23| . 

The state described by is actually quite a general 
form; it is a mixed state (unless there is equality in crite- 
rion (1401) 1 accounting for noisy processes (when /i ^ 0). 
It was proposed to provide a description of interferome- 
ters with losses [24| and the density matrix given by (j3~i|) 
represents a wide variety of states formed by non-linear 
optics experiments [23j |. 

We are concerned with inputs of the form p — 
e /(2)(9a,Pa) e ff ( 2)(>j6,Pb) j suc h f ac torizable two- mode states 

are defined by the product of the Wigner functions de- 
scribing each mode. For the output after the beamsplit- 
ter to retain factorizability requires that the output is in 
the form of the product of Wigner functions describing 
each mode, i.e. 

W a (q a ,p a )W b (q b ,p b ) ->W^q a ,p a )Wl(q b , Pb ). (41) 

-DO 

Since we are writing the Gaussian Wigner functions in 
an exponential form, factorizability is again equivalent 
to there being no cross terms in the exponent of the 
output — the Wigner function of the output is of the form 

g/(2) .P* ) gff(2) (9b ,Pb) _ 

Thus impinging a Gaussian state of the form (fM)) into 
one input port of a beam splitter and another Gaussian 
into the other mode (with corresponding parameters p! 
and t') the output is an exponential and we find that the 
coefficients of the cross terms q a qb, PaPb, qaPb are 



2rt(p' + pi* + t') 2rt(p + p* + r) 



2rt(-p' - (j,'* + t') 2rt(-p + -p* + t) 



l\2 



72 _ 4 | M , 

2irt(p' - p'*) 
r' 2 -4|^| 2 



t 2 -4|/x| 2 
2irt{p + -p*) 
t 2 -4|^| 2 ' 



(42) 
(43) 
(44) 



respectively (the coefficient of the other cross term p a q b 
is identical to that of q a p b )- Setting these equal to zero 
yields two solutions: (p = 0, r = 0) and (p = p',r = r'); 
the first is rejected as it is unphysical (it fails the positive 
semidefmite criterion Eq. (|40[> : the second is valid. Note 
that it places no restriction on the value of zq, thus when 
considering Gaussian input fields an arbitrary displace- 
ment has no effect on the factorizability of the outputs. 



IV. SUMMARY 



We discussed the role of correlated states in quantum 
state engineering, in particular in conditional measure- 
ment schemes. It was highlighted that the correlations 
required to implement such a scheme need not be entan- 
glement of the modes but rather that classical correla- 
tions suffice. This led us to consider what class of states 
input to a beam splitter exit as factorizable states; this 
class of states possess neither classical correlations nor 
quantum entanglement. By considering the action of the 
beam splitter on arbitrary inputs we found that only in 
the case that the two input states are identical Gaus- 
sians, apart from an arbitrary displacement applied to 
each mode, do we obtain an output that is factorizable. 
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